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Single Impurity Anderson Model with Coulomb Repulsion between 
Conduction Electrons on the Nearest-Neighbour Ligand Orbital 
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We study how the Kondo effect is affected by the Coulomb interaction between conduction 
electrons on the basis of a simplified model. The single impurity Anderson model is extended to 
include the Coulomb interaction on the nearest-neighbour ligand orbital. The excitation spectra 
are calculated using the numerical renormalization group method. The effective bandwidth on 
the ligand orbital, D"'^ , is defined to classify the state. This quantity decreases as the Coulomb 
interaction increases. In the D'^^ > A region, the low energy properties are described by the 
Kondo state, where A is the hybridization width. As D"'^ decreases in this region, the Kondo 
temperature Tk is enhanced, and its magnitude becomes comparable to A for D"^ ^ A. In the 
region, the local singlet state between the electrons on the / and ligand orbitals is 

formed. 

KEYWORDS: Coulomb repulsion between conduction electrons, numerical renormalization group method, mag 
netic excitation, single particle excitation, charge excitation 



Recently, much attention has been paid to the Kondo 
effect for systems with Coulomb interaction between the 
conduction electrons (CCE).^"^^ It was reported that 
a small CCE enhances the Kondo temperature Tk. A 
correlated host has usually been given by the Hubbard 
model, and the studies have been carried out by means of 
the perturbation theory on the CCE because its influence 
on the Kondo effect can be treated only by such proce- 
dures. Therefore it is not clear how the CCE affects the 
magnetic impurity problem when it is not weak. Detailed 
properties of the low energy excitation are also unknown, 
and thus it is not clear whether low energy properties can 
be described as the Kondo state. In this letter we restrict 
ourselves to a simplified model but study the effect of 
CCE in detail. We choose a model which includes the 
Coulomb interaction between conduction electrons only 
on the nearest-neighbour site of the magnetic impurity. 
By this simplification, change in the bulk host states 
caused by CCE is excluded from the consideration in 
the strict sense. But the effects through the local change 
around the impurity can be studied in detail. We expect 
that gross features caused by the CCE can be extracted 
from the studies of the present model. We calculate the 
dynamical excitation spectra such as the single particle 
and the magnetic excitation without restricting ourselves 
to the weak CCE cases. From these studies we show that 
the low energy properties are well described as the usual 
Kondo state with increased effective Tk when the effec- 
tive bandwidth of the conduction electron, D'^^ , is larger 
than the hybridization width of the /-electron, A. As 
the strength of CCE increases, D'^^ gradually decreases. 
When D'^^ becomes comparable to A, the effective Tk 
also becomes comparable to A. The excitation spectra 
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change their characteristics rapidly when the strength of 
the CCE increases further. 

In this letter we present the calculation of the dynam- 
ical excitation spectra for the following model Hamilto- 
nian. 



H = Ha + Hi^, 
Ha = H{ + Hhyh + He, 

H{ = Ef Ufa 
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Hhyh = -7^ 2^(/o^Ck<T + h.c), 
k.cr 
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(4) 
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Equation (2) represents the single impurity Anderson 
model (SIA) using the standard notation. Energies Ef 
and are measured relative to the Fermi level, E-p. 
The Coulomb interaction between electrons on the lig- 
and orbital is given by the term Hi^. The operator tilo- 
is the occupation number on the ligand orbital on the 
nearest-neighbour site, and is defined as ni,a = c^^ci^a 
with clct = (1/V^) Sk "-kcr- Here the quantities El and 
J7l are the energies of the orbital and the Coulomb in- 
teraction on it, respectively. We restrict ourselves to 
the electron-hole symmetric case, i.e., 2ef -\- Uff = 
and 2sl + = 0. Hereafter, we parameterize the 
hybridization intensity, V, by the hybridization width, 
A = TrV^/2D. In this letter, we use the numerical 
renormalization group (NRG) method to calculate the 
single particle and the magnetic excitation spectra. All 
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Fig. 1. The excitation spectra for various Ui, values in the cases of Uf{/D = 8.0 X 10~^ and A/D = S.Qtt x 10"'*. From (a) to (e), the 
abscissa is the excitation energy and ranges from 10~^ to 10* on the logarithmic scale in units of hybridization A. The value of C/l 
is 0.0 (a), 4.0 (b), 8.0 (c), 11.0 (d) and 14.0 (e). In these figures, the lines show the single particle excitation spectra of /-electron 
7rApf/Aff (solid line) and that of c-electron on the ligand orbital pi^/Nf (broken line), magnetic excitation spectra x'^{/^{ ^ 10-0 
(one-dot-dash line) and charge excitation spectra x"f/^f ^ 0-2 (two-dot-dash line). Here A^f is the degeneracy factor, A^f = 2. The 
DOS of the conduction electrons for Hi^ + model, pi^o/Nf, is also plotted as a dotted line on the same scale. Figure (f) shows x'^f 
on the linear energy scale in units T^™' . 



results are obtained using the discretization parameter 
A = l.S.^-'We consider a band with a constant density of 
states which extends from —D to D for simphcity.^-* 

Figure 1 shows the excitation spectra for various cases 
of L^L- In these figures, parameters Us and A are fixed 
at Us/D = 8.0 X 10-3 and A/D = S.Ott x lO-^.io) For 
C/l = 0, this parameter set gives the Kondo regime of 
SIA, Us > A. 

First, we examine the J7l dependence of PlO) which is 
the density of states (DOS) on the nearest-ligand orbital 
for a fictitious model setting as A = 0. It is indicated by 
the dotted lines. For A = 0, the Hamiltonian (1) is de- 
composed into two independent parts, Hf and i?L + ^^c- 
The part TJl + He corresponds to the Wolff model of the 
magnetic impurity. The quantity pLo is normalized to 
have an integrated intensity equal to 1. Therefore, the 
intensity at w = should be 0.5 according to Friedel's 
sum rule. As seen from Fig. 1(a) the present calcula- 
tion fulfills this condition with about 97% accuracy, and 
seems to be sufficient for obtaining qualitative conclu- 
sions. 

When ?7l = 0, pLo is given by the rectangular DOS 
which extends from —D to D as shown in Fig. 1(a). As 
seen from Fig. 1(b), pLo shows a three-peaked-structure 
for the case U-l/2 > D. The peak which appears 
in the energy range |w/A| < 10^ corresponds to the 
Kondo resonance of the c-electron for the i?L + He 
model. Furthermore, satellite peaks are also found at 



uj/A ^ ±2.0 X 10'^. These energies correspond to the 
atomic excitation, uj/D ~ ±[/l/2. As ?7l increases, the 
width of the Kondo resonance of c-electron becomes nar- 
row. This may be regarded as the reduction of the ef- 
fective bandwidth caused by the CCE. We define the 
effective bandwidth, , as the peak position of the 
magnetic excitation of the Hi^ + He model. This quan- 
tity almost coincides with the half-width of the Kondo 
peak of pLo- We have the relation D"^ = 0.951i:'(EE Dq) 
for Ul = 0- The quantity D"^ is shown by the double 
circles in Fig. 2 as a function of J7l- 

Next, we see how the magnetic excitation spectra of 
the /-electron, Xmt, depend on J7l for A ^ 0. In Fig. 1, 
x'^f is plotted as a one-dot-dash line. Its intensity is 
obtained by multiplying the energy of its peak position, 
r^'"^ We adopt this value, T^™ , as the characteristic 
energy of the low-energy excitation, because it is lower 
than that of other excitations such as charge fluctuation. 
We note here that T^™' has a value of about 0.6Tk for the 
conventional SIA in the Kondo regime, where Tk denotes 
the Kondo temperature determined from the magnetic 
susceptibility at T = 0, X: as Tk = l/4x- 

In Fig. 1(a), T^'^ is at w/A = 5.87 x IQ-^. As C/l in- 
creases, T^™^ initially shifts to the high-energy side (from 
Fig. (a) to (d)) and returns to the low-energy side (from 
(d) to (e)). We note that the maximum value of T^™"* 
has magnitude comparable to A ((d)). 

We replot the magnetic excitation spectra in Fig. 1(f). 
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The energy is scaled by T^™^ , and the spectra are normal- 
ized by the peak heights. The spectra almost coincide up 
to values oiUi^/D < 8. This means that the magnetic 
excitation has an almost universal shape, and therefore 
the low energy properties can be described as the Kondo 
state characterized by the energy scale As men- 

tioned later, this feature can be seen when the condition 
jjeS y ^ holds. On the other hand, the line shape of 
Xmf rapidly changes to narrow at around D'^^ ~ A when 
J/l increases further. 
The quantity, as a function of Ui,/D for var- 

10~" 



10- 



10 



-2 



H 

Q 

'a 

Q 



10- 



10- 



UfpO.O 



P 

1.6 



i n 0.0 
" » • • • • 



, HHfc A A * 

8.0x10" 



10" 



10" 




" I'exio"^ 



Fig. 2. The Kondo temperature defined as the energy of the peaJc 
position of the magnetic excitation spectra of the /-electron, 
/D, and the effective bandwidth defined as the energy of 

the peak position of the magnetic excitation spectra of the c- 
electron for the Hi^ + He model, /D, as a function of Ui^/D. 
The double circles indicate D'^*^ /D, and the other open symbols 
denote T^^/D for A/D = S.Ovr x IQ-^ cases. The solid symbols 
indicate T^^/D for A/D = S.Ott x 10"* cases. Curves are a 
guide to the eye. 



ious Us cases is shown in Fig. 2. It increases as J7l 
increases in the D^^ > A region, while it decreases 
gradually in the D^^ < A region. The magnitude of 
the maximum value is comparable to A both in the ex- 
treme Kondo regime {Us ^ A) shown by the solid dia- 
mond symbols and the valence fluctuation (VF) regime 
{Us "C A) shown by the solid circle symbols. In the 
extreme large J7l region, i.e., D*^^ <C A region, T^™^ 
decreases as T^^/D oc 16A/7r(i7L + Us). 

In the previous studies, the increase of in the 

weak CCE region is ascribed to the renormalization of 
the Kondo exchange interaction, J.^'^'^^ It is modified 
from J to J(l + 7) with the renormalization factor, 
J (xU-l/D, when one uses the perturbation theory in the 
lowest order of the Coulomb interaction. To check the 
relation, we examine the Ui, dependence of Tj^\ The 
modification of J leads to the relative change of T^^'^) 



p(m) 

[k 

p(in) 



exp 



ttUs 7 
8A 1 + 7 



(7) 



Here we use J = 8A/nUs for the symmetric SIA and 



is the value of T^""' for J7l = 0. The quantity 7 is 
obtained from eq. (7) 



7 = 



8A 

ttUs 



log 



(m) 



(m) 



KO 



- 1. 



(8) 



Using the results shown in Fig. 2, we found that the 
linear relation 7 = aUh/D holds well in the range < 
U-l/D < 0.6. Here a itself is a function of A and Us- The 
value of a is given as {Us/D, a) = (0.8, 2.6 x 10~^) and 
(1.6,3.6x10-1) for A/£) = S.Ott xlQ-^, (8.0x 10-^ 5.4x 
10-3) and (1.6x 10-2, 8.8x IQ-^) for A/D = S.Ottx IQ-"*. 

In order to examine the larger [/l region, we replot 
the results of Fig. 2 in terms of the relation T^™^/A vs 
{A/ Dp) ■ {Dq/D^^ - 1) as shown in Fig. 3. Here Dp i s the 
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Fig. 3. The Kondo temperature defined as the energy of the peaJc 
position of the magnetic excitation spectra, T^™^/A, as a func- 
tion of (A/Do)-(I?o/0''*-l) for three values of ;7fi/A. The dot- 
ted lines show the fitted results using 7' = a' {A/ Do) (Do / D"^ — 
1) and eq. (7). See the text for detail. Here a' = 3.2 {Ua/A = 
80/37r) and 6.0 (160/37r). 



value of D'^^ for U^ = 0. The quantity Dp/D^^^ corre- 
sponds to the enhancement factor of the specific heat of 
the c-electron for Hi, + model. If D'^^ is expressed as 
D'=^ = Dq exp(-7r;7L/4D), the quantity {Dq/D''^ - 1) is 
proportional to Ul/D for small Ui,/D. When the value of 
Us/ A is fixed, the data points for various A cases are al- 
most on a common curve. We try to reproduce the curves 
by using the quantity 7' = a' {A/ Do){Do/ D''^ - 1) and 
substitute it for 7 in eq. (7), where a' is the fitting param- 
eter. The enhancement factor 7 in eq. (8) is proportional 
to {A/Do){Do/D^« - 1) in {A/ Do){Do/ - 1) < 0.1 
region. We determine the value of a' within this region. 
The fitted results arc plotted by dotted lines in Fig. 3, 
and they reproduce T^'""* qualitatively even in the region 
D^^ ~ A. 

Let us return to Fig. 1. The solid and the two-dot- 
dash lines are the single particle excitation spectra (SPE) 
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of the /-electron, pf, and the charge excitation spectra, 

x"f, respectively. The intensity of pf is illustrated by 
multiplying the factor ■n/S./N^. Therefore, the intensity 
at w = should be 1.0 in Fig. 1(a) for the case of C/l = 0. 
The present calculation fulfills this condition with about 
97% accuracy. In Fig. 1(a), the Kondo resonance appears 
in the energy region |a;/A| < 10~^, and the atomic-like 
excitation appears at w/A ^ ±4.0 which correspond to 
Lo ~ ±t/fj/2. The peak position of x'li almost coincides 
with that of the atomic-like excitation. The DOS on 
the nearest-ligand orbital, pL, show a depression in the 
energy region |w/A| < 10^^ as shown by the broken 
line. It behaves as pL lo"^ in this region. Two peaks 
appear on each side of the depressed region. They have 
shoulder and peak structures which are caused by the 
Kondo resonance at oj/A ~ ±10~^ and by the atomic 
excitation at a; ~ ±f7ff/2, respectively. 

When C/l increases, the width of the Kondo resonance 
increases in the D'^^ > A region as seen from Figs. 1(b) 
and 1(c). This increase can be ascribed to the enhance- 
ment of the effective exchange coupling, and is consistent 
with the increase of T^^"^ . The intensity of pi at uj ~ 
gradually decreases when C/l increases. These facts 
indicate that the hybridization coupling between / and 
ligand orbital is enhanced by the Coulomb interaction. 
The atomic excitation peak is shifted towards higher en- 
ergy by a small amount and the spectral shape becomes 
narrow. 

In the case of 1)°* < A, the magnetic excitation has 
a rather sharp peak at the energy defined as T^'"\ The 
intensity of pf in the energy region |w| < T^™^ is very 
small as shown in Fig. 1(e). On the other hand, pf and pL 
have small peaks at about |a;| ~ The energy of the 

atomic excitation is not significantly changed. The exci- 
tation energy xj^^ is roughly given by 16A/7r(C/L -I- Us), 
as already noted. This energy corresponds to the ex- 
change coupling between the electrons on the / and 
the nearest-ligand orbital. These facts may indicate 
that the / state and its nearest-neighbour state form 
a local spin singlet pair when the charge fluctuation is 
strongly suppressed by the C/l term. This local singlet 
pair couples weakly to band states in the outer region. 
As noted previously, x'mf rapidly changes its line shape 
around the cross-over region from the Kondo singlet re- 
gion (D^^ > A region) to the local singlet one (D'^^ < A 
region). 

We note that the cross-over can also be seen from 
the analysis of the flow chart of the renormalized en- 
ergy level (FCEL) in the NRG.^^^ We have analyzed 
the FCEL of the present model, and foimd that the fixed 
point is the local Fermi liquid state in both cases, but 
the difference lies in the asymptotic area to the fixed 
point. In the > A case, the low energy levels are 
similar to that of the conventional SIA, i.e., they are 
classified as the ordinary Kondo state. However, in the 
^eff ^ ^ case, the energy levels indicate that the local 
spin singlet pair and the conduction electrons are almost 
independent of each other. 

In summary, we have examined the impurity Ander- 
son model with CCE on the nearest-ligand orbital. The 



interaction C/l reduces the effective bandwidth for the 

c-f hybridization, D'^^ . The characteristics of the sys- 
tem depend on whether the effective bandwidth is 
larger or smaller than A. In the D'^^ > A region, the 
characteristic energy of the magnetic excitation, which 
is defined as T^™"* , increases until its magnitude becomes 
comparable to A when C/l increases. The low energy 
properties are expected to be given as the Kondo state. 
The increase of T^™-* is approximately expressed as the 
enhancement of the hybridization. When C/l increases 
further into the D""^ < A region, the magnetic excita- 
tion energy decreases as 16A/7r(C/L -|- Ug), refiecting the 
formation of the spin singlet pair state from the electrons 
on the /-orbital and its nearest-neighbour orbital. 

The present model includes CCE only on the ligand 
orbital to which the /-orbital has the direct hopping 
matrix. ^"^^ However, we expect that the present results 
will be generalized to the model with CCE for the band 
states, if one concentrate on the effects due to the modifi- 
cation of the hybridization process. The quantity D'^^ of 
the present model should be interpreted as the effective 
bandwidth of such a band CCE model. The calcula- 
tion including the self-consistent procedure based on the 
d = oo model to extract a true effective bandwidth will 
be attempted in the near future.^'''' 
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